Why should primordial perturbations be in a vacuum state? 
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y^Q : Abstract 
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In order to calculate the power spectrum generated during a stage of inflation, we have to specify 

(N ' 

q ■ the quantum state of the inflaton perturbations, which is conventionally assumed to be the Bunch- 

<D . 

Davies vacuum. We argue that this choice is justified only if the interactions of cosmological 

perturbations are strong enough to drive excited states toward the vacuum. We quantify this 

efficiency by calculating the decay probabilities of excited states to leading order in the slow- 
r- • 

roll expansion in canonical single- field inflationary models. These probabilities are suppressed by a 

OO '. 

slow-roll parameter and the squared Planck mass, and enhanced by ultraviolet and infrared cut-offs. 

(N 

For natural choices of these scales decays are unlikely, and, hence, the choice of the Bunch-Davies 
■ vacuum as the state of the primordial perturbations does not appear to be warranted. 

6 
to 

in 
> 



'Electronic address: armen@phy.syr.edu 



I. INTRODUCTION 

Temperature maps of the cosmic microwave background [l| quite conclusively show that 
the universe was homogeneous on super-horizon scales already during its early stages. If 
the age of the universe at that time had been infinite or very large, this observation would 
not pose any particular problem; the universe would have had plenty of time to become 
homogeneous and isotropic by recombination. However, this happens not to be the case. 
As we follow the evolution of the universe back in time from decoupling, we soon encounter 
a time when the energy density becomes Planckian. At the Planck time, our description 
of gravity looses is validity and we are forced to impose initial conditions on the state 
of the universe. There are then two logical ways to explain the isotropy of the cosmic 
microwave background. Either the universe was already homogeneous at the Planck time, 
or homogeneity at recombination was the result of a dynamical process in the very early 
universe. 

The first explanation is possible but somewhat problematic. Effects above the Planck 
energy density could be responsible for homogenizing the universe, but since our current 
theories do not apply (or are not sufficiently developed) in that regime, those processes are, 
by definition, beyond the reach of present theoretical scrutiny. It is also possible that the 
universe was simply born in a homogeneous and isotropic state, but since this is a matter 
of initial conditions, any other initial state is possible as well. 

The second way to explain cosmic homogeneity requires a stage of inflation Q , the expla- 
nation currently favored by our present understanding of the universe. One of the successes 
of inflationary models is that they do not need particularly fine-tuned initial conditions to 
explain the properties of our universe. Say, in chaotic inflation [3| the stage of accelerated 
expansion that leads to a spatially flat homogeneous and isotropic universe is an attractor 
solution of the background equations of motion |4| . A patch of a few Hubble volumes at the 
Planck time containing a sufficiently homogeneous scalar field will start inflating almost 
independently of the initial value of the scalar field or its time derivative. In short, infla- 
tionary models are successful because they are able to dynamically explain the seemingly 
fine-tuned initial conditions of the old Big-Bang model (see jf| for different point of view.) 

Perhaps an even more important success of some inflationary models is the explanation of 
the origin and properties of primordial perturbations [7]. But alas, in order to explain these 
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properties, one has to assume that perturbations were initially in a very particular state, 
the Bunch-Davies vacuum. The Bunch-Davies vacuum is defined to be the quantum state 
devoid of quanta at asymptotic past infinity. If a perturbation mode spent an infinite time 
in the short-wavelength regime, the Bunch-Davies vacuum would be the natural state for 
the perturbations to be in, since excitations would have had an infinite time to decay into 
the ground state. This is not what happens though. Consider for instance the evolution of a 
fixed comoving perturbation mode. Because its physical momentum monotonically increases 
back in time, at some point in the past the unknown higher-dimensional operators in the 
theory that we safely neglect at low momenta actually become dominant. At that point the 
theory breaks down and we are not able to follow the evolution of the mode past that point; 



this is the well-known trans-Planckian problem jg]. 

Since we do not know the theory that describes cosmological perturbations at high mo- 
menta (early times), we have to specify the quantum state of any perturbation mode at a 
sufficiently early, but finite time. 1 The earliest possible moment at which this is meaningful 
is the time at which the physical momentum of a mode equals the cutoff of the theory. Again, 
one could simply postulate that trans-Planckian effects placed the state of the perturbations 
in the Bunch-Davies vacuum (or any other similar state), or that the latter is the natural 
state for the perturbations to be born into. But, as before, there is no justification of any 
of these choices, since they involve regimes beyond the reach of our theories. 

Therefore, the question we are going to address here is why cosmological perturbations 
should initially be in the Bunch-Davies vacuum, a state we — perhaps arbitrarily — happen to 
single out as the state of the perturbations. In order to satisfactorily answer this question, 
we have to study the predictions of our theories in the regime where they apply, namely, 
below the cut-off. Thus, we shall consider whether processes below the Planck scale are able 
to drive any initial state of the perturbations at cut-off crossing to a state that resembles the 
Bunch-Davies vacuum. In order to do so, we shall calculate transition probabilities between 
excited states of the perturbations and the Bunch-Davies vacuum. If these probabilities are 
high, these states effectively decay into the Bunch-Davies vacuum during cosmic history, and 
we are entitled to use the latter to predict the properties of primordial perturbations. In 
analogy with the homogeneity problem, we can then say that the Bunch-Davies vacuum is an 



1 For the sake of clarity we frame these heuristic arguments in the Schrodinger representation. 
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attractor for the quantum evolution of the perturbations. On the other hand, if transition 
probabilities are small, excited stays do not decay into the Bunch-Davies vacuum, which 
looses its status as the dynamically preferred state of the inflaton perturbations. Guth 
and Pi 9] made a related observation long time ago in the context of new inflation 
They noted that given the strength of the self-couplings of the inflaton needed to obtain a 
realistic spectrum of density fluctuations, thermalization would not have had time to occur 
by the beginning of inflation, so that the assumption of a thermal state for the inflaton 
perturbations might not be necessarily justified. 

II. FORMALISM 
A. Background 

Let us consider an inflating spatially flat Friedman-Robertson- Walker universe, 
ds 2 = a 2 (if) [—df] 2 + dx 2 ]. In conformal time coordinates, time runs from 77 = —00 in the 
asymptotic past to rj = in the asymptotic future. We shall assume that inflation is driven 
by a single canonical scalar field if slowly rolling down its potential. In this case, the infla- 
tionary stage can be characterized by the (dimensionless) slow-roll parameter 

where H = a' /a 2 is the Hubble "constant" and a prime denotes a derivative with respect 
to conformal time. During slow-roll inflation, the slow-roll parameter is small and nearly 
constant. We shall work to lowest non- vanishing order in the slow-roll expansion throughout 
this article. In this approximation, e is exactly constant and small. It follows then that, to 
leading order, the scale factor equals 

(2) 



as in spatially flat de Sitter. 
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B. Perturbations 



In order to describe cosmological perturbations in such an inflating universe, we choose 
a generalized comoving gauge, in which the perturbed metric reads 



ds 2 = a 2 (r]) [-(1 + 28N)drf + exp(2C)S ij (dx i + a5N i dr ] )(dx j + a5N j dr))] 



(3) 



and the scalar field perturbations vanish, 5<p = 0. The perturbation variables SN and 
SN l are constrained, and can be expressed in terms of £ and its derivatives to any desired 



order in perturbation theory, see 



10| . Note that, to linear order, ( describes the curvature 



perturbation in slices of constant scalar field [111 ] . For reasons that will become apparent 
below, we do not consider tensor perturbations. 

It is going to be convenient to expand the curvature perturbation £(77, x) in a discrete set 
of Fourier modes. We hence assume that perturbations live in a compact toroidal universe 
of side L and comoving volume V = L 3 , where they can be expressed as 
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£(77, x) = (k(v) exp(i k ■ x), with k = —j-n. 



(4) 



At the end of the calculation we shall take the continuum limit L — > 00. Note that our 
normalization implies that is also dimensionless. Because the variable ( is real, (j: = ^. 



C. Quantization 



Cosmological perturbations can be quantized by conventional canonical methods 12]. 
Substituting the perturbed metric into the Einstein-Hilbert action and keeping terms up to 
quadratic order, one obtains the free action [13] 

V 



So 



k 



v k v -k 



(5) 



where the quantization variable v? is given by 



2eaM P (- k 



(6) 



and Mp is the reduced Planck mass, Mp = (8ttG) 1 (recall that to lowest order in the 
slow-roll expansion e is constant.) In order to quantize the theory, we write the Heisenberg 
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operators as v% = v^{r])a^+v*^{r])a} It then follows that the modes v% satisfy the differential 
equation equation 

-£).* = <>. (7) 

We normalize the mode functions by imposing the condition V ■ {v^v*~ — v~Vj>) = i, which 
implies [a^, at] = 1. 

D. States 

Choosing properly normalized mode functions vArj) amounts to selecting a set of ladder 
operators a%. We can use the latter to construct a complete set of states in the Hilbert space 
of the mode k, 

fl 



l%> = ^H%>, wh ere o s |0 s > = 0. (8) 



'Til 

Following convention, we shall select here mode functions that approach the Minkowski 
space behavior vt oc e~ lkv in the asymptotic past, 



-ikrj 
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implicitly defining the Bunch-Davies vacuum. The choice of mode functions is not particu- 
larly important, as any other set also leads to a complete set of states that we can use to 
express any state in the theory. Consider for instance the set of mode functions 

Vu = A m + B^ (10) 

where \A^\ 2 — \BA 2 = 1 (this condition ensures that the mode functions are properly nor- 
malized.) Because the states \n^) associated with the mode functions v% form a complete 
set in the Hilbert space of the mode k, we can express the vacuum state |0^) associated with 
the mode functions in terms of the Bunch-Davies states \n%) [jjj], 

i fo I / B* \ " 



where we have assumed that is real. However, one cannot single out the Bunch-Davies 
vacuum by arguing that |0^) contains quanta \nt) and is hence non-empty, since it follows 



from the analogous relation 



1 V2n\ ( Br\ n w 



that the Bunch-Davies vacuum contains quanta \n%). 

This arbitrariness in the choice of vacua simply reflects that in a free th eory , specially if 



time translation symmetry is broken, there is no preferred quantum state 15]. Indeed, at 
this point there are only two natural conditions that we might impose on the state of the 
perturbations during an inflationary stage \ip). Because the background we are studying is 
homogeneous and isotropic 2 , we shall require that the quantum state of the perturbations 
itself be invariant under translations and rotations. It follows then that the expectation value 
of ( is a constant, and because ( itself is a perturbation around a homogeneous background, 
we can demand the tadpole condition (ip\(\ip) = 0. Note that these conditions are satisfied, 
for example, by any state with a definite number of excitations of each mode k, 

k 

where the occupation numbers only depend on the magnitude of the wave vector, = n k . 
E. Interactions 

Interactions between the different modes of cosmological perturbations arise from their 
self-couplings, and from the couplings to additional matter fields that reheating requires. 
The former are model-independent, in the sense that couplings can be determined uniquely 
for given a set of slow-roll parameters, whereas the latter are model-dependent, in the sense 
that couplings to matter are hardly constrained by the requirement of a successful reheating. 

The cubic self-couplings of cosmological perturbations have been calculated in [101 ]. ex- 
tended to fourth order in jig ], and generalized to k-fields in [itJ (see also [3]). To lowest 
order in the slow-roll expansion, they are determined by the action 



5S = M 2 P VJ2 J dr ia 2 e 2 

ki,k 2 



Ck 1 Ck 2 Ck 3 ~ ' hCk 1 Ck 2 Ck 3 - 2 2 Cfci Cfc 2 Cfc 3 



^2 



(14) 



where k 3 is fixed by momentum conservation, k 3 = —k\ — k 2 . Because of the unconventional 
normalization of the variable (, one can not directly read off the cubic coupling constants 
from the action. Noting that the Heisenberg operators are related to the "canonically" 



2 Isotropy is broken by the periodic boundary conditions we are imposing on the spatial sections of the 
universe. This breaking can be ignored on small scales. 
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normalized by equation (jBJ), we find that interactions are suppressed by y/e/Mp, as ex- 
pected from gravitational couplings. Interactions between scalars and tensors are suppressed 
at least by an additional factor of e, and can thus be neglected to leading order in the slow- 
roll approximation. Quartic couplings are suppressed by an additional factor of the Planck 
mass and are hence further suppressed at momenta below the Planck scale. Note by the way 
that interactions vanish in the de Sitter limit e — > 0. In the theory at hand, the interaction 
Hamiltonian is 5H = —SS. 

The couplings in equation (I14p satisfy criteria that guarantee that expectation values 
of cosmological perturbations are dominated by the contributions of sub-horizon modes 
20, l2l(]. The two terms that contain time derivatives are "safe", whereas the term with two 
spatial derivatives is "dangerous." For our purposes it is going to be convenient to use safe 
interactions only. Using the free equation of motion ([7]), integrating by parts, relabeling 
the fields and invoking momentum conservation, we arrive at the equivalent interaction 3 (to 
leading oder in the slow roll expansion) 

5S = M P VJ2 Jd V a*e*.(l + 2^( kl ('J' kr (15) 

Note that the non-derivative inflaton couplings to scalar fields usually postulated to be 
responsible for reheating do not satisfy Weinberg's criteria of convergence 2^] , although 



it can be shown [21] that they do not lead to large quantum corrections to cosmological 
correlation functions. 



F. Expectation Values 

Inflationary predictions about the properties of primordial perturbations typically involve 
the expectation value of an appropriate observable. The observable most widely considered 
is the power spectrum V(k), which is implicitly defined by the equation 

m\k,o)ak',om = ^6 u ,v{k). (ie) 

3 In the interaction picture, operators satisfy the free equations of motion. The actions are equivalent if 
the boundary terms obtained upon integration by parts vanish, which is actually not the case at the 
asymptotic future rj = 0. Because decay probabilities should not depend on the interactions long after a 
mode has left the horizon, we shall ignore this boundary term. 
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The power spectrum only depends on a single comoving momentum. For simplicity we shall 
consider here observables 0^(77, fc) of this type, where the argument k labels the operator 
and the subscript H identifies the observable as an operator in the Heisenberg picture. 

As an example, let us consider the power spectrum for an excited state (fT51) in a free 
theory. Using equations fl§D, ©, (TTO]) and rotational symmetry we find 

l + 2n k H 2 j . 

PW = (17) 

Clearly, in a free theory the power spectrum depends on the state we choose, a dependence 



that underlies most of the "trans-Planckian" effects discussed in the literature 



22] . 



G. Interaction Picture 

In order to calculate the expectation value of an observable Oh(t], k) in the presence of 
interactions, it is convenient to work in the interaction picture. Let us split the Hamiltonian 
of cosmological perturbations into free and interaction parts, TC(rj) = Ti-o(r]) + 5H(r]). The 
Hamiltonian of cosmological perturbations is explicitly time dependent. In the interaction 
picture, operators carry the free time evolution, and the interacting Hamiltonian (in the 
interaction picture) generates the time evolution of state vectors. It follows then that 

&\0 H (ri, kM) = (iP\U}{r), —T)Oi(r), kp^rj, -T)|^>, (18) 

where — T is an arbitrary time in the past, at which the interaction picture is introduced, 
and 

Ojfa, k) = U \ V , -T)O h (-T, k)U ( V , -T), 
U (n,Vo) =Pexp [-ifl T dri"Ho(v')], Uj( V , Vo ) = Pexp [ - i $\ drfSH^rf)] . (19) 

Typically, we are interested in evaluating the expectation value of an observable like 
the power spectrum long after the mode has left the horizon during inflation. It is hence 
appropriate to consider the asymptotic future 77 = as the upper limit in the integrals. 
The choice of the lower limit T in the integrals is more ambiguous. Conventionally, one 
considers the asymptotic past (T — > 00), but there are reasons to believe that this limit is 
not physically realistic. Indeed, several arguments suggest that general relativity is just an 
effective field theory, valid up to an ultraviolet cut-off A uv , above which the theory looses its 



validity. Because the physical momentum of the mode k monotonically increases into the 
past, the earliest time T k at which our theory can faithfully describe it is 

1 A uv =► T k = ^, (20) 



a{-T k ) uv " Hk 

where the subindex indicates that this time depends on the comoving momentum k. Since 
— Tfc is also the time at which we introduce the interaction picture (where Heisenberg and 
Schrodinger representations agree), the state of this mode in the Heisenberg picture can be 
regarded as the initial state of the mode in the Schrodinger representation. 

In spite of the previous objections, one can also formally consider the limit T — ► oo. 
In this limit, Ui(rj, —T) does not converge, and one is forced to introduce a regularization 
factor. Usually this is achieved by adiabatically switching off the interactions in the past, 
that is, by replacing the coupling constants e by ee £r> and letting e — > at the end of the 
calculation. Whereas this procedure is justifiable in scattering experiments, where the wave 
packets representing the different particles do not overlap in the past, it is less clear what it 
means in an inflationary spacetime. This formal prescription is useful nevertheless, because 



a theorem due to Gell-Mann and Low 23] guarantees that (at least in Minkowski space) 
the limit selects the interacting eigenstates of the theory. Hence, we shall also consider the 
class of states defined by the relation 

m = uf\-T,-ooM, (21) 

where the superscript (e) denotes that the couplings constant e of the mode k has been 
replaced by e exp(e^), with e = 1/T k . Because the operator Jj\ is unitary, the new set 
of states also form a complete set in the Fock space of the different modes. The factor 
exp(er]) in the coupling constant turns off the interactions at times r\ < — so at these 
early times the quantum state of each mode is an eigenstate of the free theory. In that 
respect, equation (12111 describes states that at time 77 = — T k are given by \n k ). By slowly 
turning on the interactions however, the prescription (T2~Tj) places the system closer to the 
actual eigenstates of the full interacting theory, whatever they are. Note that equation (12 ip 
is merely a formal definition of a class of states, and does rely on any particular assumptions 
about the interacting theory above the cut-off A uv . 

To conclude this section, let us express the expectation value of a Heisenberg picture 
operator in a slightly different form. By inserting two expansions of the identity into equation 

10 



(TT8j) we find that 

^|^(r/)|^) = ^|(0|^(^-r)|^)| 2 (0|O / (^)|0), (22) 

where we have assumed that the non-diagonal matrix elements of the interaction picture 
observable vanish. This is in fact the case for observables like the power spectrum, for which 
the expectation value of the interaction picture operator is diagonal and only depends on the 
state of the mode k. In summary, there are two ways to calculate the expectation value of 
an operator Ouij]-, k). One can calculate it directly using equation (1181) . or one can evaluate 
the transition amplitudes {<p\Ui(rj, — T)\ip) and substitute into equation (|2~2"|) . It turns out 
that the evaluation of the transition amplitudes will help us to understand the dependence 
of the expectation value on the state of the perturbations. Note that the expectation values 
of our approximate eigenstates of the interacting theory, equation fl2T|) . are given by 

(^|OH(r/)|^) = El^l f/ f £) ^- °)l^| 2 ^l^^)l^' ( 23 ) 



where \<f>), \ip) are the "free" states of the form (|T3"|) . 

III. TRANSITION PROBABILITIES 

To zeroth order in perturbation theory, the expectation value of the operator O is given by 
{iP\Oh(v^)\4') ~ (i J \^i( r l: Because interaction picture operators carry the free time 

evolution, this is just the result one would get in the free theory, which is the approximation 
commonly used in calculations of the primordial spectrum of cosmological perturbations. 
From equation fl22l) . this approximation is justified as long as the forward transition proba- 
bility 

P=MU I (v,-TM)\ 2 = \T\ 2 (24) 

is close to one. In order to calculate P it is going to be convenient to employ the unitarity 
of the operator U[. Since Ujlli = I, it follows that 

P = l-P tot , where P tot = ^ MU^r], -T)|^)| 2 (25) 

is the total decay probability of the state \ip). Therefore, just as in field theories in Minkowski 
space, all we have to do to calculate the forward scattering probability P is compute the 
decay probability P tot . 
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We shall calculate in the following P tot to first order in perturbation theory. In order to 

— * 

do so, we shall assume that all the modes but k are in the Bunch-Davies vacuum, 

w = \ n k) n i°*> ( 2g ) 

To further simplify and reduce the amount of terms we have to consider, we shall restrict 
our attention, without significant loss of generality, to the subset of final states 

, AV^+l)- 1 alAAW (2?) 
,t „t 



n k al 3 al 2 a k \i/j) : 

where k 2 is an arbitrary momentum and, because of momentum conservation, /c 3 = k — k 2 . In 

— * 

a transition from \ip) to the number of quanta in the mode k increases by one, so these 
processes represent mode excitations. In a transition from \ip) to the number of quanta 

— * 

in the mode k decreases by one, so these processes represent emission or decay. For notational 

— * — * 

convenience, we shall identify in the following k = k±, which we shall interchangeably utilize. 



A. Decay into 0_ 

To lowest order in the slow roll expansion (and perturbation theory), the decay amplitude 
into the state is given by a single integral over time of the appropriate products of mode 
functions. Evaluating the integral, and defining the "energy" change AA; = k 2 + k 3 — ki we 
arrive at the transition amplitude 



AV^MpV k x k 2 k 3 



1 - e~ iAkT ■ (1 + iAkT) l-e~ iAkT 
Kl {Kk/W + AJfc/Jfcx 



(28) 



where the dimensionless coefficients K\ and K 2 are given by 

( 1 1 1W , , k z k\ k\k 2 \ , T ^ k 2 k 3 k 2 k 3 

We are interested here in calculating transition amplitudes to leading order in the short- 
wavelength parameter k{T . In the limit k{T ^> 1, the expressions inside the square brackets 
that explicitly depend on Ak approach functions that enforce energy conservation. The 
amplitude of processes for which Ak/ki > (A^T)" 1 is suppressed at least by a factor of 
(kiT)^ 1 , and the amplitude of processes for which 

AA; 1 
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is enhanced by the inverse of the same quantity. Because our effective theory description is 
valid up to h\T = k\T k = A uv /H, the latter are the transitions that yield the leading order 
in the ultraviolet cut-off. For energy conserving transitions, equation (I30p . the transition 
amplitude is 

1 e^ 2 H \~n^ 



— ifi • {k{T k f + 3iK 2 - k x T k 



(31) 



4 V 1 I 2 M P V kxhh 

It thus appears, from equation ([201 . that the leading term in the transition amplitude is 
proportional to A^ v . However, this is misleading, since for small values of Ak/ki, 



— and K ~ ( fc i ~ k ^ + 
h 2 (h - k 2 )k\k 2 



K^K- and K 2 M ^- ^IT^i . (32) 



As a result, because of equation (1301) . we find that the contribution proportional to K\ is 
at most of the same order as the one proportional to K 2 ; the transition amplitude is hence 
proportional to A uv . 

Armed with this knowledge, we can readily estimate what the decay probability into a 
given state \<f>-) approximately is. Rather than expressing it in terms of comoving momenta, 
it is going to be instructive to express it in terms of the physical momenta involved in the 
decay. For energy conserving processes with k\/k 2 of order one, it is given by 



Ufc 1 1 ek\ 



64 a?Vk x k 2 k z M 



(33) 

phys 



The label "phys" denotes that the momenta in the equation are the physical momenta, 
a ■ fcphys = k. Note that the factor l/(a 3 V) is the physical volume of our compact universe 
at the same time. What is surprising about equation fl33|) is that the Hubble constant has 
dropped out of the probability; the only dimensionful quantity that enters the equation is 
the Planck scale. Hence, there is no suppression of the decay for modes outside the Hubble 
radius, or for decays into those modes, and the only place where the Hubble parameter 
plays a role is the allowed energy change A/c p hys < H. If all the momenta are of the same 
order k phys w A uv , the decay probability is proportional to e • (A uv /Mp) 2 . Note that the 
dimensionless factor Vkik 2 k% w Vk\ is roughly the number of states with momentum less 
than k\. 

In order to calculate the total decay probability, we shall need to precisely count the 
number of states with energy k 2 and azimuth 9 (the z axis of our coordinate system points 
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along the direction of k = ki.) In the continuum limit, the density of states in an infinitesimal 
interval around those values is 

P= dk 2 -d(cos6) = 4^ 2 ' {3 ) 

It is hence painfully obvious from equations (134]) and (|28|) that the probability density 
p ■ IT 2 ) is infrared divergent. Indeed, as k 2 jk\ — > (which implies k% ^ k\ — k2Cos6), we 
find i^i — > 1 — cos# and - ► k\jk 2 . Because in the same limit Ak/k — ► (1 — cos8)k 2 /ki 
approaches zero, to leading order in an infrared expansion (and to leading order in short- 
wavelengths) we get 4 

2 , 9e^fc H 2 k\ 2 /or^ 

P ' |T -' ^U^Mlh' Tk ' (35) 
We encounter a similar divergence in the limit k% — > 0. Since transition probabilities are 

symmetric under the interchange k 2 ^ k^, we can restrict our attention to the limit of small 

&2 (by "small" we mean k 2 <C fci.) 

To regularize the total decay probability we introduce, on top of the already discussed 

ultraviolet scale, an infrared cut-off at comoving momenta k- m f. Because there are infrared 

divergences in both limits k 2 — > and k 2 — > fei = fc, the regularized total decay probability 

is given by 

k-kinf 

P tot = y dk 2 d(cos6) p ■ Pi^f_. (36) 

We shall calculate the leading logarithmically divergent contribution to the total decay 
probability. In doing so, it is important to remember that the approximation ( |35l) is valid 
as long as equation (1311 holds. Hence, integrating the decay probability up to k 2 = 1/T k we 
finally arrive at 



P t 



9rik eA 



2 

uv 



"" 16jr 2 Ml 



l0g lr-F7l +o(1) 



(37) 



where the subleading terms of order one arise from the remaining portion of the integral, 
and the subdominant terms in the infrared expansion that we have dropped. We have also 
multiplied the resulting decay probability by two, to account for the infrared divergence 
at small k 3 . Notice that in the derivation of equation f[3T|) we assumed that > k^, 



4 Because of the infrared divergence, one has to be careful in approximating 
here. 



l-cxp(-tAfcT) 2 
(Afe) 



by 27rT(5(Afc) 
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which implies that the logarithm is positive. If the logarithm happens to be negative for 
a particular choice of k, the term of 0(1) gives the dominant contribution to the decay 
probability. 

Using equation (1251) we can calculate transition probabilities that involve the interacting 
states defined in equation (I2TT) as well. We simply substitute Ak — ie for Ak, and let T — > oo 
(taking the limit T — > oo corresponds to letting the initial time in the interacting propagator 
approach infinity.) The resulting amplitude is a function of Ak and e, where e = 1/Tjt is 
finite and depends on the time at which interactions are "turned on" . For energy conserving 
processes, equation ( 1301) . the transition amplitude is again given by equation ( 13~T1) . up to the 
factor of 1/2 multiplying K\. Hence, had we used these interacting states in the calculation, 
we would have obtained the same result. 

Infrared divergences often occur in theories that contain massless particles. It is hence not 
surprising to meet these divergences, since the mode functions of cosmological perturbations 
behave as those of massless fields. In QED for instance, the leading correction to the cross 
section for electron scattering with no soft photon carrying away an energy larger than E 
diverges 5 as E — » [24]. Because photon detectors cannot detect photons of arbitrarily low 
energies, the measurable cross section is always finite, and determined by the precision of the 
detectors involved in the measurement. Similarly, we do not expect to be able to determine 
corrections to the power spectrum that involve modes of arbitrarily long wavelengths. At 
the very least, our observations are constrained by our present horizon. Hence, identifying 
the infrared cut-off with our present comoving horizon ko, we might set 

few « k . (38) 

Equations fl37|) and fl38l) imply that the decay probability are too small to justify the assump- 
tion that perturbations are initially in the Bunch-Davies vacuum. Current measurements 



suggest that the slow roll parameter e is of the order 10~ 2 25j, whereas the validity of our 
approach requires the ultraviolet cut-off A uv not to exceed the Planck mass Mp. Hence, 
for our present horizon, we expect a decay probability P t ~ t < 10~ 3 n^. The window of 
observationally accessible modes encompasses at most five decades in /c-space, so for all 
observable modes the decay probability remains small. Let us stress that this conclusion 



5 In some cases, one can sum these divergences to all orders in perturbation theory and obtain a finite result 
even in the limit E — > 0. Such a summation however is beyond the scope of our paper. 
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crucially depends on the choice of the ultraviolet cut-off. Other choices lead to much higher 
probabilities. 



B. Excitation to <j) + 

The calculation of the transition amplitude 7^,^ + is analogous to the calculation of 
71. The result has the form of equation (|2"gj) . where, as opposed to the previous case, the 



amplitude is proportional to y/n k + 1 and Ak equals k± + k2 + k%. The coefficients K\ and 
K 2 are given by 

( 1 1 1W, , k z k\ k\k 2 \ T ^ k 2 k 3 k 2 k 3 

Again, the total transition probability is infrared divergent. In order to isolate the leading 
contributions, we consider the limit of small k 2 , for which we can approximate Ak/k ~ 2. 
Because there are no energy-conserving transitions in this case, the term proportional to 
K\ — > k\jk 2 dominates the transition amplitude. The time dependent factor multiplying K\ 
is —iT/Ak, regardless of the value of k 2 . Following the procedure described in the previous 
section we then arrive at 

n k + leA 2 u 



P+ ~ 

1 tot 



uv 







log 


( —) 


\hnij 



(40) 



64tt 2 Ml 

Note that this time, the logarithm contains the ultraviolet cut-off. 

Equations (1371) and (HOT) suggest that the vacuum (ng = 0) is as stable or unstable as its 
excited states (n^ ^ ). What is happening here is that we have actually misidentified the 
vacuum [10 ]. Let us now calculate the total decay probability of the interacting vacuum in 
equation (pTj) and see whether we recover again fT40l . To calculate amplitudes that involve 
the interacting states, we replace Ak by Ak — is and let T —>■ oo. The result is 

1 e l / 2 H jn k + l 



AV l / 2 Mp\l hk 2 k : 



(41) 



where K 1 and K 2 are given by equation (|39l) . Although the explicit dependence on T k 
appears to be gone, one should bear in mind that the factors Ak now contain e = 1/T k . 
In this case this dependence is irrelevant, since the total probability is dominated by the 
contributions from modes at small k 2 , for which AA; = k± + k 2 + k 3 — is ~ 2k\. Therefore, 



16 



integrating over k% we find 



2 



P* 



49 (n k + l)eH 
tot ~ 256tt 2 Ml 



(42) 



which shows that excitations from the interacting vacuum are suppressed with respect to 
decays from excited states by a factor of (H/A uv ) 2 . Inserting equation (j4"2l into fl22|) . we 
recover the magnitude of the one-loop quantum corrections to the power spectrum derived 



in Q. 



C. Non-standard vacua 

Because of equation (TTTT) . the results of the previous subsections can be applied to non- 
standard choices of vacua. In particular, since the state \0k) contains Bunch-Davies quanta 
| rife), our analysis of transition probabilities can be indirectly used to estimate the survival 
probability of any state that contains |0fc). Instead of following this route, it proves 
handier to compute the transition probabilities directly. 

Let us assume that the state of the mode k is |0fc), and let us calculate what is the 
probability of transition to a state \k, ^2,^3) = ^!^! 2 ^! 3 |0fc)- If the mode functions only 
contain positive frequency components, the amplitude of a transition |0fc) — ► &2, £3) is 
suppressed, because the time integrals of the mode functions lead to factors that enforce 
energy conservation, as in equation fl28l) . However, if the mode functions contain negative 
frequency components, some of these energy conservation factors contain "wrong" signs, 

J vC kl C k2 'C k3 'd V oc A* kl At 2 A* k3 ■ ( + ^ + (43) 

/ 1 - e^ kT l-e-^(l + tAfcT) x 
+ B kl A k2 A k3 { - -^== j + • • ■ , 

where Ak = k\ + ki + k 3 , Ak = —k\ + k 2 + /c 3 , and the dots denote many terms that 
we have not written down. Whereas Ak is what we would conventionally identify as the 
energy change, the structure of Ak suggests that the mode k\ has negative energy —k\. As 
a consequence, transitions from the vacuum to excited states are kinematically allowed, and 
their amplitude is proportional to the negative frequency components B^. These results 
simply show that the vacuum states |0fe) are unstable. Since the vacuum has a component 
with negative energy, it can "decay" into positive energy quanta. Note however that these 
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instabilities arise when interactions are turned on. In the absence of interactions, any state 
is stable, and for appropriately weak interactions, any state is sufficiently long lived. 

The total transition probability from the vacuum |0fc) into any state with three quanta 
can be calculated along the same lines as before. Once again, the probability is infrared 
divergent. The contribution from the modes around k 2 ~ is given, to leading order in 
A uv /Ai n f, by ^ 

P ^4^Mf l0g (^) ■ \ A ^ A ^B kl - A kl B kl B k2 \ 2 , (44) 
where we have used the interacting states defined in equation (I2ip . Note that most of the 
non-standard vacuum choices discussed in the literature satisfy A « 1 and B < H/A uv 22]. 
Hence, for these states the transition probability is at most 

P^^oJ±), (45, 



4tt 2 M 2 p ta VMnf, 

which is of the same magnitude as the probability of excitation into the Bunch-Davies state 
in equation (1421) . as the reader might have expected from equation (HI]). Given that 
constraints on the energy scale of inflation imply that H < 10~ 4 Mp [261 ]. the transition 
probability is very small; non-standard vacuum choices of this kind appear to be legitimate. 



IV. FORWARD SCATTERING PROBABILITIES 



In our previous calculations we have taken the comoving cut-offs of the theory to corre- 
spond to the physical cut-offs at the initial time rj = —T k . Because for a fixed physical scale 
A the corresponding comoving momentum k = aA increases with the expansion of the uni- 
verse, we might have underestimated the phase space available in the different transitions. 
In this section we calculate directly the forward scattering probabilities, a check that shall 
allow us to verify our estimates. The reader not interested in technical details might skip 
this rather technical section, which just confirms our previous results. 



A. Lowest Order 



The forward scattering amplitude (ip\Ui\ip) receives contributions from the Feynman 
diagram in Figure [TJ The state \nz) decays into two quanta \k 2 , k%) at time r\i and the latter 
annihilate into |nr) at time r] 2 > T]i. The amplitude for the process contains an integral 
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FIG. 1: Lowest order contribution to the forward scattering amplitude. One of the excitations 
with momentum k in the state |nr) decays at time rj\ into two modes with momenta k 2 and ks, 
which at time r] 2 annihilate into a state with momentum k = k\. 



that runs over all possible values of k 2 at time 771, which, discarding vacuum fluctuation 
diagrams, is given by 

m k 

(V>IW> ~ 1 - J dr) 2 J dru J dk 2 dcosOp ■ M*(r/ 2 )M(r/ 1 ) e - iAfc (,?2 - r ' l) , (46) 

-T -T k ini 

where, as before, Ak = k 2 + k 3 — ki, and 



M{r]) 



80(77) M P V Vhk 2 k 



1 + 2-i ] (1 — — — ) ■ (—k 2 k 3 ) ± permutations 



(47) 



For convenience, we have made the time ordering in the perturbative expansion explicit, 
because the phase space available for the decay depends on the comoving cut-off at the 
decay time rji, k — a(?7i)A uv for instance. 

For any fixed value of 771, the integral over momenta in equation (146!) is infrared divergent. 
Let us calculate again the leading contribution in an infrared expansion. For small values 
of k 2 we find 



P M* ( V2 )M( Vl 



3n k eH 2 2 
64tt 2 M 2 p 1 



+ i(l - cos 9) ■ (772 - 7/1) 



+ 



(48) 
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Therefore, substituting the leading term into equation (|46j) we arrive at 

o m k 

mM) ~ 1 " I dr]2 J dr]1 I dk 2 dc ™ e Y 2 e ~ lAHm ~ m) - (49) 

P -T -T k ini 

Since k 2 is assumed to be small, let us approximate Ak « k 2 (l — cos#) in the exponen- 
tial. Integrating the latter over the cosine gives a factor that effectively cuts off the integral 
at k 2 ~ (i] 2 — Tjx)" 1 . Therefore, the integral over momenta yields a time-dependent fac- 
tor 21og[(?72 — ^i) _1 /^mf] ? which, since (r] 2 — cannot exceed k, can be replaced by 
— log[/;; in f • (i] 2 — f]i) + k in {/k]. Integrating this expression over rji and r] 2 we finally arrive at 

ww-'-sf^K)' (50) 

in agreement with equation ( 1371) . 



B. Next to Lowest Order 



Because the leading order in the infrared expansion depends on the logarithm of a ratio 
of comoving scales, the forward scattering amplitude is insensitive to the time at which 
these comoving scales are evaluated. In order to check whether this property is preserved 
at subleading orders, we shall calculate the next to leading order in the infrared expansion. 
Substituting the subleading term in equation fHHj) into f|46l) and integrating over cos 6* we 
find 

r)2 k 

A//lJ7l/\~ 3n ^ H2 L f, ( \ >„2 [ l-e- 2 ^-^[l + 2ik 2 (r l2 -r ]l )} 

— T -T 

(51) 

Again, the rational term effectively cuts off the momentum integral at values 
k 2 > l/(rj 2 — r)i). Hence, the momentum integral yields 1/(772 ~ /71 ) , and we finally obtain, 
after integrating over the time variables 

« (52) 

also in good agreement with the subleading term in equation (J37J. 

In summary, the direct calculation of the forward transition amplitudes confirms our 
previous transition probability estimates. In a field theory in Minkowski space this agreement 
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is just a reflection of unitarity, but in our context this connection might have failed because 
comoving cut-offs are time-dependent. In the indirect approach, one integrates first over 
time, and later over comoving momenta. Since comoving cut-offs are time-dependent there 
is in principle an ambiguity in the time at which the integral over momenta has to be cut 
off. In the direct approach, one integrates over momenta at the time when transitions occur, 
which makes the comoving cut-off unambiguous. Given the quite different nature of both 
calculations, it is reassuring that both approaches agree. 

V. CONCLUSIONS 

The choice of the Bunch-Davies vacuum as the quantum state of the inflaton perturbations 
can be justified only if interactions efficiently drive excitations toward that state. In this 
paper we have calculated the decay probability of excited states to leading order in slow-roll 
approximation and to leading order in an ultraviolet and infrared expansion. The decay 
probability of a scalar perturbation mode k with nk quanta, P t ~ t , is approximately given 
by equation (157)) . Because the self-interactions of the perturbations are proportional to the 
slow-roll parameter e, the decay probability vanishes in the limit of de Sitter inflation. As a 
consequence, the closer the spectral index is to one, the more unlikely are excited states to 
decay. 

Since cosmological perturbations behave as massless fields, and because their interactions 
are gravitational, the transition probabilities are infrared and ultraviolet divergent. Inspec- 
tion of the decay probability reveals that, for natural choices of the corresponding cut-offs, 
an excited state is not likely to decay during cosmic evolution, P t ~ t < 10~ 3 for modes close 
to our present horizon. Although excited states are significantly less likely to survive than 
the (interacting) Bunch-Davies vacuum, equation (1421) . the small value of the decay prob- 
ability suggests that there is no reason to assume that primordial scalar perturbations are 
in a Bunch-Davies-like vacuum state. Since the infrared divergence is only logarithmic, this 
conclusion is rather insensitive to the infrared cut-off. On the other hand, because the prob- 
ability quadratically diverges with the ultraviolet cut-off, our conclusions crucially depend 
on the latter. If we extrapolate our results to trans-Planckian momenta, A uv > Mp / y^e, we 
obtain decay probabilities of order one. 

Although we have not considered tensor modes here, similar results should apply for 
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gravitational waves. Because their interactions are suppressed by additional powers of the 
slow-roll parameter e, excited states are even more unlikely to decay, and non- vacuum sig- 
natures in the primordial spectrum of tensors should be even more prominent. 

We should also emphasize that the estimate in equation flTTI) is a lower bound on the 
decay probability. We have considered only the model-independent interactions that nec- 
essarily arise in any slow-roll inflationary model. Model-dependent couplings to the decay 
products of the inflaton should boost the decay probabilities, but given our lack of con- 
straints on the strength or nature of these couplings we cannot take them properly into 
account without additional assumptions. It is also important to keep in mind that our 
results only apply to canonical, single-field inflationary models. For non-canonical scalar 
fields, there are additional model-independent interactions that might lead to significantly 
higher decay probabilities. This is actually what we exp ect a p riori, since k-inflationary 



models [27|, |28|, |29| lead to stronger non-Gaussianities |17|, [18|, |30[ |. 

The ability of inflationary models to make definite predictions about the origin of struc- 
ture hinges on a dynamical justification for the choice of the quantum state of the pertur- 
bations. Our calculations suggest that the very same feature that guarantees the success 
of inflationary models, the existence of a weakly coupled slow-roll regime, renders the pre- 
dictions of the theory sensitive to the arbitrary choice of the initial quantum state of the 
perturbations. 
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